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Algorithm for muon electromagnetic shower reconstruction 

S. Mangano'' on behalf of the ANTARES coUaboration 
'^NIKHEF, Kruislaan 409, 1098 SJ Amsterdam, The Netherlands 

The ANTARES neutrino telescope is presently being built in the Mediterranean Sea at a depth of 2500 m. The 
primary aim of the experiment is the detection of high energy cosmic muon neutrinos, which are identified by the 
muons that are produced in charged current interactions. These muons are detected by measuring the Cerenkov 
light which they emit traversing the detector. Sometimes a high momentum muon produces electromagnetic 
showers. The subject of this paper is a method to reconstruct these showers which includes several steps: an 
algorithm for the fit of the muon track parameters, preselection of detected photons belonging to a shower, and 
a final fit with the preselected detected photons to calculate the electromagnetic shower position. Finally a 
comparison between data obtained with that part of the detector that is currently in operation and simulations 
is presented. 
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1. Introduction 

The ANTARES [1] detection principle relies 
on the observation of Cerenkov (CK) light emit- 
ted by relativistic charged particles in water. 
The emitted CK photons are detected by a 3- 
dimensional grid of 900 photomultiplier tubes 
(PMTs) arranged in 12 detection lines at a depth 
of 2500 m in the Mediterranean Sea. The detec- 
tor will be completed in early 2008. Since the 
end of January 2007, the detector has 5 operat- 
ing detection lines. Amongst others first data of 
atmospheric downward going muons have been 
measured, which are up to 500 m long tracks 
with energies larger than 100 GeV. This gives the 
unique opportunity to study the average number 
of showers per track length for high energy muons 
passing through sea water. 

The muon track is reconstructed from the ar- 
rival times of the CK light detected by the PMTs, 
whose positions are known. The muon energy 
loss mechanism in water arc ionization, e+e~ pair 
production, bremsstrahlung, and photonuclear 
interactions. Above several hundred GeV the 
muon energy loss is dominated by the latter three 
effects [2], which are discrete and are character- 
ized by large energy fluctuations. Furthermore, 
the average muon energy loss due to these effects 
is proportional to the energy of the muon. 



In the following a method to reconstruct dis- 
crete electromagnetic (EM) showers is presented, 
which gives additional information about a muon 
track (e.g. the shower multiplicity) as compared 
to existing algorithms only measuring the direc- 
tion and position of the muon track. 

2. Identification of muon EM showers 

The challenge of EM shower reconstruction is 
to distinguish CK photons from EM shower pho- 
tons. The EM shower light has two key charac- 
teristics: it is produced in one point on the muon 
path and it arrives, in general, delayed compared 
to the CK light. 

The first characteristic distinguishing CK light 
from EM shower light is illustrated in Fig. [TJ 
The CK photons are emitted continuously along 
the muon path, whereas EM showers are emit- 
ted stochastically and discretely. Furthermore, 
the number of photons from EM showers is pro- 
portional to the deposited energy. By project- 
ing the 3-dimcnsional position and arrival time 
of all detected photons onto the muon axis, the 
identification of EM showers is reduced to a one- 
dimensional problem. This transformed informa- 
tion is presented by a one-dimensional histogram, 
in which the EM shower position appears as a 
peak. So, the peak indicates where on the muon 
trajectory the EM shower is produced. 
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Figure 1. A muon emits Cerenkov photons con- 
tinuously on its trajectory in water with some- 
times a discrete electromagnetic shower (upper 
panel). The emission points of the photons from 
the muon trajectory are calculated from the mea- 
sured 3-dimensional position and arrival time in 
the PMTs. These emission points are presented 
by a one- dimensional histogram with its Z-axis 
aligned with the muon axis. In the histogram the 
peaks identify the electromagnetic shower position 
(lower panel). 



sidered relative to the expected arrival time of a 
direct CK photon from the reconstructed muon 
track. This calculation can be simplified by a 
rotation of the coordinate system. The M hits 
x'^ = {x'^,y'i,zl,t'^ with i = 1,....,M of an event 
are rotated by the rotation matrix: 



Ri9, 



where 9 and (j) are given by the muon direction. 
The situation after the rotation Xi — R{9, 0)a;^ 
of the i-th measured hit a^ = (xi^yi, Zi,ti) as 
well as the rotated muon space-time coordinates 
(x, y, z, t) is presented in Fig. ^ 
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3. Algorithm for shower reconstruction 

The algorithm consist of three steps using dig- 
itized PMT signals which are referred to as hit^. 
The first step selects muon candidates. Then, 
hits are preselected which originate from one dis- 
tinct shower. The third step is a x^ -minimization 
yielding the space-time position of the shower 
within the preselected hits. 

3.1. Selection of shower hits 

This section presents the mathematical frame- 
work which is schematically illustrated by Fig. [1] 

Muon events are selected by using an existing 
muon reconstruction algorithm described in [3], 
which provides an estimate of the direction and 
position of the muon at some fixed time. For 
all hits in one muon event, the hit time is con- 



^ A hit occurs when a PMT signal crosses a threshold 
corresponding to 0.3 of a single photon-electron amplitude. 



Figure 2. Schematic view of the detection of the 
Cerenkov light and the randomly emitted shower 
light. The thick line represents the muon, the thin 
line the path of the Cerenkov light and the thin 
dotted line the path of the shower light. The muon 
goes through point {x,y,z,t). The Cerenkov light 
is emitted at an angle 9ck with respect to the 
muon track at point Zf^ and is detected by a 
PMT as a hit in point (xi,yi, Zi,ti). The ran- 
domly emitted shower light is emitted at point Zi 
and is detected by the same PMT, but later in 
time. 



The expected CK arrival time (tf^) of the pho- 
ton is calculated by (see Fig. [2|): 



("" = ( + 



c V 



tan 9cK ^ c sin 9ck 



.(1) 



Roma International Conference on Astro-Particle physics 



where c is the speed of hght in vacuum, n is 
the refraction index of water, Ock is the CK 
angle and r^ is the perpendicular distance be- 
tween the muon trajectory and the PMT given 

by r^ = y^{xi - a;)2 + (y^ - y)2. 

3.1.1. Calculation of the photon emission 
point from the muon 

In order to suppress random background hits, 
a selection is made of hits for which the differ- 
ence between the calculated arrival time tf'^ of 
CK photons and the measured arrival time is 
in the interval —t,nin < ti — tf^ < t^ax- Further- 
more, this interval is subdivided in two unequal 
intervals, because photons from EM showers ar- 
rive in general later than CK photons. The early 
time interval contains CK photons and is given 



by \t^ - 
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The late time interval 
■ix contains mostly EM 



shower photons. The value for tmin is given by 
the timing resolution of the PMT and the disper- 
sion of light in water, whereas tmax is defined by a 
value where the number of signal hits approaches 
the background level. 

The emission location of photons Zf^ from the 
early time interval is calculated under the CK an- 
gle assumption and is given by (see Fig. [2|): 

tan 6cK 
For the late time interval the unknown photon 
emission location Zj is given by (see Fig. [2]): 

Zj — z 
t^ = t+ — + 



+ -\/{xi 
c 



xf + {y,~yY + {z^-Z,f,{?,) 



where the second term is the time for a muon 
moving at speed of light to reach the point Z^ 
where the light is emitted, while the third term 
is the time needed for the light in water to travel 
from Zi to the PMT. Eq. [3] corresponds to a 
quadratic equation with two distinct solutionqj 
represented by Z+ and Z^ . 

3.1.2. Peak finder 

All calculated Zf ^, Zf and Z~ positions are 
presented by a one-dimensional histogram. The 



^The late time interval ensures that two real solutions ex- 
ist. 
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Figure 3. The upper panel represents all recon- 
structed emission points of the photons emitted 
by an atmospheric muon. They are calculated 
from the measured hit information after a full 
simulation of the detector response to an atmo- 
spheric muon event with backgrounds noise. The 
middle panel shows the hits which are selected 
from the algorithm, after the peak finder step (see 
text). Only statistically relevant peaks have been 
selected. Furthermore, hits produced from each 
generated shower are represented by a different 
hatching, emphasizing that the purity of each peak 
is high. The lower panel displays the three gener- 
ated shower positions on the muon axes. 



upper panel of Fig. [3] shows an example of a full 
simulation of the detector response to an atmo- 
spheric muon event with background noise. 

These Zf^, Z^ and Z~ positions are used by a 
one-dimensional peak finder algorithm. It selects 
peaks with amplitudes higher than a threshold 
identifying only statistically relevant peaks. The 
peak finder takes into account that if the two solu- 
tion Zj and Z[^ are used in different peaks, the 
solution in the smaller peak is eliminated. The 
result of the peak finder is shown in the middle 
panel of Fig. [31 In addition hits produced from 



Muon electromagnetic shower reconstruction 



a particular shower are represented by a different 
hatching. The non-hatched hits are either pro- 
duced by CK hght or background noise. 

In the lower panel of Fig. [3]the generated posi- 
tion of the three showers are also presented. The 
agreement between the middle and lower panel 
shows that the selected peaks in this event have 
a high hit-purity and the algorithm succeeds in 
properly reconstructing the shower position on 
the muon axes. 

3.2. Reconstruction of the shower position 

The last step of the algorithm is to calculate 
the EM shower space-time position under the as- 
sumption of a point-like shower source for each 
peak. For this, a first estimate of the shower 
space-time position xl = {xs,ys, Zs,ts) as ex- 
plained in [3] is used. Then for a peak with N 
hits x'j = (a;' , y', z', i') a x^ with a fixed a can 
be introduced: 



^ 1 



i;--t.- 



^{xr - X,)2 + iy'^ - y,)2 + {^'^ - ,^)2^ (4) 



If necessary the x^ is minimized by removing 
hits until among all possible permutations of the 
hits used to calculate the shower position a satis- 
factory x^ normalized to the number of freedom 
is found. 

The result of the algorithm is that for each 
muon candidate the shower multiplicity can be 
calculated. For each shower candidate different 
quantities can be evaluated: e.g. the position of 
the shower, the closest approach between shower 
position and muon track and the number of hits 
per shower. 
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Figure 4. Distribution of the number of recon- 
structed atmospheric muon shower candidates for 
both data (points) and simulation (dotted line). 



Fig. m The agreement between data and simu- 
lation is satisfactory. For the used set of selec- 
tion cuts, the one shower muon events constitute 
about 5% of all muon events. 

In summary a method to reconstruct EM show- 
ers from a muon has been developed and applied 
for the first time in ANTARES. The essential el- 
ement of the algorithm is that the selection of 
shower hits is reduced to a one-dimensional prob- 
lem. Improvements of the algorithm as well as a 
more extensive application to data are subject of 
the future work. 



4. Conclusion 

The performance of the reconstruction algo- 
rithm for EM showers induced by muons has been 
validated in a sample of simulated multiple at- 
mospheric muon events with background noise. 
Next, the algorithm was applied to data obtained 
with ANTARES. The results of the shower mul- 
tiplicity for down going muon events is shown in 
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